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The electronic structure and magnetic properties of iron-base superconductors including lithium- 
iron-arsenic (LiFeAs) and iron-selenium (FeSe) have been studied using density functional theory 
and hybrid functional. The computed anti-ferromagnetic exchange interactions for the first (second) 
nearest neighbouring iron atoms in LiFeAs and FeSe are approximately 16 (81) meV and 6 (13) meV 
respectively. These results are in agreement with previous theoretical calculations and experiments. 
The total energies for LiFeAs and FeSe unit cell of different spin configurations on iron atoms includ- 
ing non-magnetic, spin-1, and spin-2 are also calculated. The lowest state has a spin-2 configuration 
for LiFeAs and FeSe. These calculations show the co-existence of the complicated band structure 
around the valence band and conduction band and d-orbital magnetism. The theoretical results 
presented here are particularly useful for the experimentalists working on the electronic structure 
and magnetism of LiFeAs and FeSe. 
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I. INTRODUCTION 

Cuprate-based superconductors (CuSCs) have the 
highest superconducting transition temperature achieved 
to date in HgBaCaCuO (T c - 200 K) [1]. By con- 
trast a rather new class of high-T c superconductors 
(HTSC), iron-based superconductors (FeSCs) have at- 
tacted much attention recently following the discovery 
of the high-transition-temperature superconductivity in 
LaFeAs(0,F) Q. FeSCs have the second highest transi- 
tion temperature, e.g., the recently achieved T c ~ 55 K 
in SmFeAsOo.85 [3]. It was surprising that iron atoms 
can play a crucial role in high-T c superconductors. The 
FeSCs that have received extensive experimental and the- 
oretical studies include ReFeAsO (Re=La or rare earth) 
(1 Ill-type) @, H, AFe 2 As 2 (A=Ba, Sr, Ca, etc) (122- 
type) 0, LiFeAs (FigUJi) or NaFeAs (Ill-type) [5], and 
FeSe (11-type) @ (FigQJ)). FeSCs are more tunable than 
CuSCs as one can control transition temperature not only 
by doping like in CuSC but also by applying pressure on 
parent stoichiometric compounds. 

The common structural properties among these FeSCs 
is that iron atoms form two-dimensional layers, in which 
spins on iron atoms could interact with each other 
via a super-exchange pathway mediated by non-metallic 
species [7[. This could be naturally modelled by using 
a two-dimensional Heisenberg model [7[. This feature is 
also shared by CuSC whereas copper atoms carry spins. 
However, the superconducting gap structures of FeSCs 
and CuSCs may be different [13]. However, recently 
there is intense debate on the existence of magnetism, 
especially in Li(Na)FeAs (8l4l2j. These debates are con- 
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centrated on the electronic structure of FeSCs near the 
Fermi surface. 

The most important topic in HTSC is electron-pairing 
mechanism. More generally speaking, what mechanism 
supports superconductivity to survive at such a high 
temperature? Bardeen- Cooper- Schrieffer (BCS) theory 
[lH stated that the electron-pairing mechanism in con- 
ventional superconductors is phonon-assisted thanks to 
the Fermi-liquid behaviour of electrons in metals mak- 
ing the attraction overcome repulsion near the Fermi 
surface. However, in HTSC it is broadly believed 
that phonon alone cannot lead to HTSC. Neverthe- 
less, this belief has never been proven correct carefully. 
Based on this belief, people would rather like to at- 
tribute the dominant pairing mechanism to the mag- 
netism induced by the anti-ferromagnetic (AFM) interac- 
tion [15| . The co-existence of the first-nearest-neighbour 
and second-nearest-neighbour AFM exchange interaction 
in FeSCs leads to a competition between Neel ordering 
and collinear ordering 0, EH, [I?}- Complicated band 
structure near the Fermi surface and magnetic proper- 
ties in FeSC and CuSC are crucial for the understanding 
of the mechanism of HTSC. However, most of the first- 
principle calculations performed so far are in the level of 
local density approximation (LDA) and generalized gra- 
dient approximation (GGA) and a satisfactory descrip- 
tion of the band structure taking into account electron 
correlation properly is still rare 

[ill The hybrid 

functional such as B3LYP [U and PBEO @ can provide 
a better description for electron correlation and partially 
eliminates the self-interaction error by mixing exact ex- 
change. For example, the band gap in various compounds 
can be predicted more accurately by the density func- 
tional theory and hybrid functional than those computed 
by LDA [13 • 

In this paper I report ab. initio calculations for LiFeAs 
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FIG. 1: (Colour on line.) The structures of a 2a x 2b x 2c 
super-cell of LiFeAs (a) and FeSe (b) are shown. In (a) small 
brown balls are used to label lithium atoms, large brown for 
arsenic, and small yellow for iron. In (b) large green balls are 
used to label selenium, and small yellow for iron. 



and FeSe by using density functional theory (DFT) and 
hybrid functional PBEO [25]. The rest of this paper is 
organized as the following: in ^TT] I will introduce the 
computational details, in §1111 1 will show my calculation 
results and discuss them, and in §IVI I will draw some 
conclusions. 



II. COMPUTATIONAL DETAILS 

The calculations for electronic structures of LiFeAs and 
FeSe are carried out by using DFT and hybrid functional 
PBEO [25| as implemented in the CRYSTAL 09 code j2fij|. 



The basis sets specially designed for lithium[27], iron[28], 
arsenic [29j, selenium [30[ atoms in solid-state compounds 
are used throughout all the calculations. The Monkhorst- 
Pack samplings (3l| of reciprocal space are carried out 
choosing a grid of shrinking factor to be 8 x 8 x 4 in 
order to be consistent with the ratios among reciprocal 
lattice parameters in LiFeAs and FeSe (their experimen- 
tal tetragonal structures [HI, [33| are tabulated in Ta- 
ble [Tj) . The truncation of the Coulomb and exchange se- 
ries in direct space is controlled by setting the Gaussian 
overlap tolerance criteria to 10 -6 ,10 -6 ,10 -6 ,10 -6 , and 
10- 12 0. The self-consistent field (SCF) procedure is 
converged to a tolerance of 10 -6 a.u. per unit cell (p.u.c). 
To accelerate convergence of the SCF process, all calcu- 
lations have been performed adopting a linear mixing of 
Fock matrices by 30% . 

Electronic exchange and correlation are described us- 
ing the PBEO hybrid functional [25| which is free of em- 
pirical parameters. The advantages of PBEO include a 
partial elimination of the self-interaction error and bal- 
ancing the tendencies to delocalize and localize wave- 
functions by mixing a quarter of Fock exchange with 
that from a generalized gradient approximation (GGA) 
exchange functional [25j. 

The broken-symmetry method [34| is used to localize 
opposite electron spins on each atom in order to describe 
the AFM state. The exchange couplings in Heisenberg 
model [35] is defined here as, 

H = Jinn ^ Si • Sj + J2NN ^ Si'Sj, (1) 

ijelNN ij£2NN 

and determined by 

Jinn = (^fm^ ~ ^afm)/( 85 ^ ( 2 ) 
J2NN = (^fm^ _ ^afm)/(4# 2 ) - 2Jiatat. (3) 

where ^afm anc ^ ^fm^ are defined as the total energies 
of a unit cell, in which the first nearest neighbouring spins 
are in AFM configuration or ferromagnetic configuration 
(FM), respectively. Jinn is defined as the first-nearest- 
neighbouring exchange interaction. The exchange inter- 
action with the second nearest neighbour J2atat, is cal- 
culated in a 2a x b x c super-cell. C^afm) are 
defined as the total energies with the second-nearest- 
neighbouring spins aligned (anti-aligned). These formu- 
lae are obtained based on the two-dimensional lattice 
formed by iron atoms as shown in Fig|2j 

DFT total-energy calculations have some intrinsic dis- 
advantages for computing exchange interactions, e.g., a 
comparison between two very large numbers to get a 
small exchange splitting in an all-electron local basis set 
formalism. However, the performance of the hybrid func- 
tional, e.g., B3LYP or PBEO as implemented in CRYS- 
TAL has previously been shown to provide an accurate 
description of the electronic structure and magnet ic prop - 
erties for both inorganic and organic compounds |36l-[40j . 

The spin states of iron atoms in LiFeAs and FeSe are 
still under active debate |8l4l2|. In this paper, the to- 
tal energies of the unit cell with non-magnetic, spin-1, 
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Symmetry 


Lattice constant (A) 


Atomic positions (in fraction) 


LiFeAs 


P4/nmm (129) 


a=b=3.791; c=6.364 


Li (0.25, 0.25, 0.8289); Fe (0.75, 0.25, 0.5000); As (0.25, 0.25, 0.2897) 


FeSe 


P4/nmm (129) 


a=b=3.7685;c=5.5194 


Fe(0.25, -0.25, 0.0); Se(0.25, 0.25, 0.2877) 



TABLE I: The structural parameters for LiFeAs and FeSe. 




a 



FIG. 2: (Colour on line.) A sketch for the two-dimensional 
iron layer in LiFeAs and FeSe is shown. Iron atoms are la- 
belled by blue balls. Jinn is the first-nearest neighbour ex- 
change interaction, and J2NN is the second-nearest neighbour 
exchange interaction. I also show the unit cell and super-cell 
(doubling the unit cell along lattice vector a) for the computa- 
tion of Jinn and J2NN respectively. The spin configurations 
in the unit cell and super-cell as shown here are used to calcu- 
late the total energies of AFM configuration £"afm an d ^afm 
respectively. 

and spin-2 state for iron atoms are calculated for LiFeAs 
and FeSe carefully to determine the most stable magnetic 
state. 



III. RESULTS AND DISCUSSIONS 

The calculations of a unit cell in LiFeAs and FeSe with 
4 unpaired electrons on iron atoms (spin-2 state) are cho- 
sen to illustrate the electronic structure and magnetic 
properties. 

A. LiFeAs 

I first present the calculation results for electronic 
structure and magnetic properties of LiFeAs. In total 21 
bands are plotted for both AFM and FM configurations. 
11 occupied bands are plotted for spin- up and spin-down 
in AFM configuration. 14(6) occupied bands are plotted 
for spin-up (spin-down) in FM configuration. The zero 
energy is chosen to be in the middle of the band gap at 
T point. 

In Figj3](a), (c), and (e), the band structure, density of 
states (DOS), and spin densities are presented when the 
spins on iron atoms in a unit cell are anti-aligned (AFM 



configuration). The valence band, the conduction band, 
and another two immediate higher bands for the spin- 
up and spin-down in AFM configurations and the spin- 
down in FM configuration are highlighted in red. In the 
AFM band structure it can be seen that at T point there 
is small gap ~ 1.0 eV between the valence band (VB) 
and the conduction band (CB). There is also a crossing 
between VB and CB between T and X point. The CB 
crosses two immediate higher bands as well. These band 
crossings make AFM configuration look like conducting 
as shown in DOS (Fig|3t) The largest gap between the 
VB and CB (~ 3 eV) can be found at L point. 

In Fig|3](b), (d), and (f), the band structure, density of 
states (DOS), and spin densities are presented when the 
spins on iron atoms in a unit cell are aligned (FM con- 
figuration). In the FM band structure, there is a rather 
large gap for spin- up, ranging from ~ 3.5 eV at T point to 
~ 8.0 eV at L point. By contrast, for spin-down, again 
we can see there is band crossing around VB and CB. 
For the spin-down of FM configuration (Figj3j}), the VB 
crosses CB at X point. And also there are several band 
degeneracies between T and L points. The analysis of 
the bands near Fermi surface suggests that the hybridiza- 
tion between d-orbitals of iron and 4sp-orbitals of arsenic 
atoms is dominant. Again the band crossings make FM 
DOS look like a half metal; there is clear band gap in 
the spin-up whereas no gap in spin-down. The Mulliken 
charge on lithium is ~ +0.7|e|, on iron ~ +0.2 |e|, and 
arsenic ~ — 0.9|e| for AFM and FM configurations. The 
Mulliken spin densities on irons with spin anti-aligned 
are ~ ±3.1/i£ for AFM configuration. The discrepancy 
between this value and the expect one, i.e., 4/i# may be 
due to the hybridization between iron and arsenic. For 
FM configuration, the Mulliken spin densities on irons 
are ~ 3.4/i^ and on arsenics are ~ 0.6/iB. This indi- 
cates there is a strong hybridization between iron and 
arsenic. The Mulliken populations for d-orbitals on the 
iron atom are also calculated and tabulated in Table ITT1 
The populations of d z 2, d yzi and d x i_ y i are slightly less 
than one because the hybridization with Asp- orbital on 
arsenic reduces the populations a bit. The populations of 
d xz and d xy are almost 0.5 and this implies one unpaired 
electron may be in a super position between d xz and d xy . 
The relatively large population of 4sp-orbitals on arsenic 
suggests a stronger hybridization in FM configuration. 

The projected densities of states (PDOS) of LiFeAs in- 
clude the projections to the 2sp-orbitals of lithium atoms 
(red), the 4sp-orbitals of iron atoms (green), and d- (blue) 
and Asp (violet)-orbitals of arsenic atoms. These show 
that lithium seems to be almost irrelevant and <i-orbitals 
on iron contribute much to the bands near VB and CB. 
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LiFeAs Mulliken populations 


Orbital 


d z 2 


dxz 


dy Z 


d x 2 —y 2 


dxy 


Px 


Py 


Pz 


FM Populations 


0.89 


0.45 


0.70 


0.79 


0.52 


0.23 


0.20 


0.10 


AFM Populations 


0.81 


0.11 


0.74 


0.75 


0.40 


0.02 


0.11 


0.08 



TABLE II: The Mulliken populations for d-orbital of iron atoms and 4sp-orbital of arsenic atoms are shown. 
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FIG. 3: (Colour on line.) The band structure (a,b), density of states (c,d), and spin densities (e,f) of LiFeAs in AFM and 
FM configurations respectively are shown. As seen in the electronic structure, AFM configuration looks conducting, whereas 
FM looks like half- metal because of the band crossing around VB and CB. The projected density of states is multiplied by a 
factor of 10 to show the orbital contribution, where 2sp of lithium is in red, Asp of iron in green, 3d of iron in blue, Asp of As 
in violet, and the total one in black. In spin densities, spin-up is in red, and spin-down in green. 
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From Fig|3] (e) and (f), it can be seen that the spin 
densities spread out to arsenic atoms when spins on iron 
atoms are aligned, and this implies the importance of the 
hybridization between iron and arsenic. The hybridiza- 
tion to the non-magnetic atoms may be able to shed some 
new light on the construction of an effective model for the 
explanation of superconductivity in FeSCs. 

The total energies per unit cell for the non-magnetic 
and spin-1 states are also calculated and both of them 
are higher than that of the spin-2 state by 6.51 eV and 
2.61 eV respectively. So the spin-2 state for iron atom 
is the most stable state. The exchange interaction in 
LiFeAs Jinn is ~ 16 meVs and J2NN is ~ 81 me Vs. The 
result that J2NN is larger than Jinn might be due to 
different bond angles for the first nearest neighbour and 
second nearest neighbour. The prediction of an AFM 
ground state in this work is qualitatively in agreement 
with recent theoretical results in Ref.[l8[. In Ref.[18[, it 
was also concluded that there is an AFM ground state 
in LiFeAs. However, the details of the band structure 
in Ref.[18j, which largely depends on the functional, are 
different from the work presented here. The previous 
work in Ref.[19[ suggests that the absence of AFM or- 
dering and points to a triplet pairing mechanism. Based 
on my results, the system will prefer a collinear order- 
ing because J2NN > Jiatat/2, i.e., each spin has two 
ferromagnetically coupling neighbours and another two 
anti-ferromagnetically coupled ones. Therefore There is 
some similarity between the work in Ref.[l9| and my 
work. However, my work disagrees with Ref.[22| be- 
cause they predicted that AFM configuration is higher 
in energy than non-magnetic one. This might be due 
to the exclusion of the exact exchange in the functional 
the authors used in [22j. My work is consistent with the 
experiments that have shown that there may be strong 
anti- ferromagnetic fluctuations [ll|, [HI . 



B. FeSe 

The band structure, DOS, and spin densities for AFM 
(FM) configuration of FeSe are shown in FigHk, c, and 
e (Fig Ob, d, and f) respectively. In total 21 bands are 
plotted for both AFM and FM configurations. 11 occu- 
pied bands are plotted for spin-up and spin-down in AFM 
configuration. 14(6) occupied bands are plotted for spin- 
up (spin-down) in FM configuration. The zero energy is 
chosen to be in the middle of the band gap at T point. 

From the band structures of FeSe (FigHJ we can see 
that both AFM and FM configurations are insulating al- 
though they have different band gap. This is qualitatively 
different from that in LiFeAs. Once again, the d-orbitals 
of iron atoms contribute a lot to the bands near the Fermi 
energy. From the spin densities, it can also been seen that 
the hybridization between iron and selenium in FM con- 
figuration. The Mulliken charge on iron is ~ +0.3 \e\ and 
on selenium ~ — 0.3 |e| for both AFM and FM configura- 
tions. The Mulliken spin densities for irons are ±3.6/i# 



in the AFM configuration. This value is close to the 
expected, i.e., 4/i#. For FM configuration, the Mulliken 
spin densities for irons are 3.7/i# and for selenium 0.3/ig. 
This also indicates the hybridization between iron and 
selenium. The VB and VB-1 for the spin-down in FM 
configuration are also highlighted in red to illustrate the 
band degeneracy. The populations of <i-orbitals on iron 
and 4sp-orbitals on selenium are tabulated in Table Hill 
The populations of d z 2, d xzi d yzi and d x 2_ y 2 are slightly 
less than one suggesting these orbitals are singly occupied 
although the hybridization reduces the populations a bit. 
We can also see the hybridization between d-orbitals on 
iron and 4sp-orbitals on selenium in FM configuration. 

I also compare the total energies per unit cell of FeSe 
for different spin configurations including spin-0, spin-1, 
and spin-2. The total energies of spin-0 and spin-1 are 
higher than that of spin-2 by 5.80 eVs and 1.92 eVs, re- 
spectively. This indicates that the spin-2 state is also 
favoured by FeSe in this tetragonal structure. The ex- 
change interaction in FeSe Jinn is approximately 6 meVs 
and J2NN is approximately 13 me Vs. This is qualitatively 
in agreement with the results for FeSe in Ref.[17j. 



IV. CONCLUSION 

DFT calculations combined with hybrid functional 
PBEO has been performed for LiFeAs and FeSe tetrahe- 
dral structure. The band structures, densities of states, 
and spin densities are shown and discussed. From these 
calculations, the AFM configuration of LiFeAs there is 
band crossing between the valence band and conduction 
band. By contrast, the AFM configuration of FeSe is in- 
sulating. And the FM configuration of LiFeAs show a 
half-metallic behaviour while FM band structure in FeSe 
is insulating. The fundamental difference between the 
band structures of LiFeAs and FeSe might be able to 
shed some light on the comparison between the exper- 
iments on them. As we have seen, the co-existence of 
localized magnetic moments and the complicated band 
structure around VB and CB might be the key to answer 
the question whether LiFeAs is non-magnetic. 

The total energies of different spin configurations in- 
cluding non-magnetic, spin-1, spin-2 are also carefully 
compared for LiFeAs and FeSe, which shows that both 
of them favour spin-2 state. The exchange interactions of 
the first nearest and second nearest neighbours are also 
calculated, and the computed values qualitatively agree 
with recently experiments and theoretical results leading 
to a two-dimensional collinear ordering rather than Neel 
ordering because J2NN > Jl * N . We can also see that 
from these calculations that the complicated band struc- 
tures around VB and CB co-exist with collinear ordering 
induced by the AFM exchange interactions between local 
magnetic moments on iron atoms. 

These calculations are particularly useful for the ex- 
perimentalists who are concerned with the electron trans- 
port, magnetic properties, and electronic structure of the 
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FeSe Mulliken populations 


Orbital 


d z 2 


dxz 


dyz 


d x 2 —y 2 


dxy 


Px 


Py 


Pz 


FM Populations 


0.89 


0.86 


0.86 


0.90 


0.07 


0.10 


0.10 


0.10 


AFM Populations 


0.90 


0.81 


0.82 


0.89 


0.05 


-0.02 


0.02 


0.0 



TABLE III: The Mulliken populations for d-orbital of iron atoms and 4sp-orbital of arsenic atoms are shown. 




FIG. 4: (Colour on line.) The band structure (a,b), density of states (c,d), and spin densities (e,f) of FeSe AFM (left) and FM 
(right) configurations respectively are shown. Notice that the band structure suggests that AFM configuration is insulating 
with a band gap about 2.5 eV by contrast to LiFeAs, whereas FM configuration is half-metal. PDOS are multiplied by a factor 
of 10 to show the orbital contribution, where Asp of selenium is in red, 3d of iron in blue, Asp of iron in green, and the total 
one in black. In spin densities, spin-up is in red, and spin-down in green. 
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Ill-type and 11-type FeSCs. In the future, I would like tronic structure information provided here to explain su- 
to build a rather novel effective model based on the elec- perconductivity in FeSCs. 
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